Extragalactic Sources Research Articles

The relationship between the turnover frequency and photo-ionization in radio sources

ABSTRACT We investigate the connection between the turnover frequency in the radio spectrum, νTO, and the rate of ionizing ultraviolet photons, $Q_{\rm {H{\small I}}}$, in extragalactic sources. From a large, optically selected sample, we find νTO to be correlated with $Q_{\rm {H\,{\small I}}}$ in sources that exhibit a turnover. The significance of the correlation decreases when we include power-law radio sources as limits, by assuming that the turnover frequency occurs below the lowest value observed. However, the power-law fit sources are less well sampled across the band and so these may just be contributing noise to the data. Given that the observed νTO–$Q_{\rm {H{\small I}}}$ correlation is purely empirical, we use the ionizing photon rate to obtain the electron density in a free–free absorption model. For each of the constant, exponential, constant plus exponential (Milky Way), and spherical models of the gas distribution, there is also an increase in the turnover frequency with ionizing photon rate. Furthermore, for a given gas mass, we find that the turnover frequency is anticorrelated with the scale-factor of the gas density. While other mechanisms, such as ageing electrons or synchrotron self-absorption, may be required to reproduce the spectral indices, for an exponential scale-factor similar to the linear size, this simple free–free absorption model reproduces the turnover-size correlation seen in radio sources.

Cosmology-independent Photon Mass Limits from Localized Fast Radio Bursts by using Artificial Neural Networks

Abstract A hypothetical photon mass, $m_{\gamma}$, can produce a frequency-dependent vacuum dispersion of light, which leads to an additional time delay between photons with different frequencies when they propagate through a fixed distance. The dispersion measure--redshift measurements of fast radio bursts (FRBs) have been widely used to constrain the rest mass of the photon. However, all current studies analyzed the effect of the frequency-dependent dispersion for massive photons in the standard $\Lambda$CDM cosmological context. In order to alleviate the circularity problem induced by the presumption of a specific cosmological model based on the fundamental postulate of the masslessness of photons, here we employ a new model-independent smoothing technique, Artificial Neural Network (ANN), to reconstruct the Hubble parameter $H(z)$ function from 34 cosmic-chronometer measurements. By combining observations of 32 well-localized FRBs and the $H(z)$ function reconstructed by ANN,we obtain an upper limit of $m_{\gamma} \le 3.5 \times 10^{-51}\;\rm{kg}$, or equivalently $m_{\gamma} \le 2.0 \times 10^{-15}\;\rm{eV/c^2}$ ($m_{\gamma} \le 6.5 \times 10^{-51}\;\rm{kg}$, or equivalently $m_{\gamma} \le 3.6 \times 10^{-15}\;\rm{eV/c^2}$) at the $1\sigma$ ($2\sigma$) confidence level. This is the first cosmology-independent photon mass limit derived from extragalactic sources.

How Well is the International Celestial Reference System Maintained in Official IAU Implementations?

The International Celestial Reference Frame (ICRF) based on the VLBI-derived positions of 608 extragalactic radio sources was adopted by the International Astronomical Union (IAU) in 1998 as the first realization of the International Celestial Reference System (ICRS). Later, in 2009 and 2020, two extended ICRF versions, ICRF2 and ICRF3, respectively, were released. The latter is adopted by the IAU as the current implementation of the ICRF in the radio band. In the meantime, the Gaia mission delivered three versions of the optical ICRS realization in 2016, 2018, and 2022 with an accuracy similar to that achieved by VLBI-based ICRF. The Gaia-CRF catalogs were linked to the ICRF under no-rotation conditions and thus may suffer from ICRF systematic instability if the latter is substantial. In this work, a new analysis was performed to assess the long-term stability of radio and optical ICRS realizations. Based on the 16-parameter vector spherical harmonics expansion of the differences between the three ICRF catalogs, it can be concluded that the mutual orientation between them is at a level of a few tens of microarcseconds, while the components of the glide vector and E 2,0 term are several times greater. A comparison of the three Gaia-CRF catalogs with the ICRF3-SX showed that for the latest Gaia-CRF catalog, Gaia-CRF3, all rotational and deformation components are below 20 μas except for the E 2,0 term, which is several times greater. For both ICRF and Gaia-CRF catalogs, the evolution of the source position errors is also tracked.

Fast computation of the observer motion effects induced on monopole frequency spectra for tabulated functions

Context. Various methods have been studied to compute the boosting effects produced by the observer peculiar motion, which modifies and transfers the isotropic monopole frequency spectrum of the cosmic background radiation to higher multipoles. Explicit analytical solutions for the spherical harmonic expansion coefficients were already presented and applied to different types of background spectrum, strongly alleviating the computational effort needed for accurate theoretical predictions. The frequency spectra at higher multipoles are inherently led by higher-order derivatives of the monopole spectrum. Provided that it can be well described by analytic or semi-analytic functions, the computation of its transfer is not affected by numerical instabilities when evaluated at the required level of numerical accuracy. Instead, monopole frequency spectra described by tabulated functions are computed with a relatively poor frequency resolution in comparison with the Doppler shift, which necessitates interpolation of the tabular representation. The spectra are also affected by uncertainties related to intrinsic inaccuracies in the modelling or in the related observational data as well as to limited accuracy in their numerical computation. These uncertainties propagate and increase with the derivative order, possibly preventing the trustworthy computation of the transfer to higher multipoles and of the observed monopole. Aims. We study methods to filter the original function or its derivatives and the multipole spectra, to mitigate numerical instabilities, and to derive reliable predictions of the harmonic coefficients for different cosmic background models. Methods. From the analytical solutions, and assuming that the monopole spectrum can be expanded in Taylor’s series, we derive explicit expressions for the harmonic coefficients up to the multipole ℓmax = 6 in terms of monopole spectrum derivatives. We then consider different low-pass filters: prefiltering in Fourier space of the tabular representation; filtering in both real and Fourier space of the numerical derivatives; interpolation approaches; and a dedicated method based on amplification and deamplification of the boosted signal. We study the quality of these methods when applied to suitable analytical approximations of the tabulated functions, possibly polluted with simulated noise. These methods are then applied to the tabulations. Results. We consider two very different types of monopole spectra superimposed to the cosmic microwave background: the (smooth) extragalactic source microwave background signal from radio-loud active galactic nuclei and the (feature-rich) redshifted 21 cm line, and present our results in terms of spherical harmonic coefficients. The direct prediction of these coefficients can be noisy at ℓ > 1 or, depending on the uncertainty level, even at ℓ ≤ 1. Without assuming a functional form for the extragalactic background spectrum, the Gaussian prefiltering coupled to the sequential real-space filtering of derivatives allows us to derive accurate predictions up to ℓ ∼ 6, while a log–log polynomial representation, which is appropriate over several decades, gives accurate solutions at any ℓ. Instead, it is difficult to characterise the 21 cm line model variety, and so it is relevant to work without assumptions about the underlying function. Typically, the prefiltering provides accurate predictions up to ℓ ≃ 3 or 4, while the further sequential filtering of the derivatives or the boosting amplification and deamplification method improves the results up to ℓ = 4, while also allowing reasonable estimations of the spectrum at higher ℓ. Conclusions. The proposed methods can significantly extend the range of realistic cosmic background models manageable with a fast computation, beyond the cases characterised a priori by analytical or semi-analytical functions. These methods require only an affordable increase in computation time compared to the direct calculation via simple interpolation.

Merger-driven Growth of Intermediate-mass Black Holes: Constraints from Hubble Space Telescope Imaging of Hyper-luminous X-Ray Sources

Hyper-luminous X-ray sources (HLXs) are extragalactic off-nuclear X-ray sources with luminosities exceeding the theoretical limit for accretion onto stellar-mass compact objects. Many HLXs may represent intermediate-mass black holes (IMBHs) deposited in galaxy halos through mergers, and the properties of the stellar cores surrounding HLXs provide powerful constraints on this scenario. Therefore, we have systematically built the largest sample of HLX candidates with archival Hubble Space Telescope (HST) imaging (24) for the first uniform population study of HLX stellar cores down to low masses. Based on their host galaxy redshifts, at least 21 (88%) have stellar core masses ≥ 107 M ⊙ and hence are consistent with accretion onto massive black holes from external galaxies. In 50% of the sample, the HST imaging reveals features connecting the HLXs with their host galaxies, strongly suggesting against the background/foreground contaminant possibility in these cases. Assuming a mass scaling relation for active galactic nuclei and accounting for an estimated contamination fraction of 29%, up to ∼60% of our sample may be associated with IMBHs. Similar to previously known HLXs, the X-ray luminosities are systematically elevated relative to their stellar core masses, possibly from merger-driven accretion rate enhancements. The least massive stellar cores are preferentially found at larger nuclear offsets and are more likely to remain wandering in their host galaxy halos. The HLX galaxy occupation fraction is ∼ 10−2 and has a strong inverse mass dependence. Up to three of the HLX candidates (12%) are potentially consistent with formation within globular clusters or with exceptionally luminous X-ray binaries.

Aliasing from Galactic Plane Setting in Widefield Radio Interferometry

Measurements with widefield radio interferometers often include the near-infinite gradient between the sky and the horizon. This causes aliasing inherent to the measurement itself and is purely a consequence of the Fourier basis. For this reason, the horizon is often attenuated by the instrumental beam down to levels deemed inconsequential. However, this effect is enhanced via our own Galactic plane as it sets over the course of a night. We show all-sky simulations of the Galactic plane setting in a low-frequency radio interferometer in detail for the first time. We then apply these simulations to the Murchison Widefield Array to show that a beam attenuation of 0.1% is not sufficient in some precision science cases. We determine that the noise statistics of the residual data image are drastically more Gaussian with aliasing removal, and explore consequences in simulation for cataloging of extragalactic sources and 21 cm Epoch of Reionization detection via the power spectrum.

Constraints on the origin of the radio synchrotron background via angular correlations

Abstract The origin of the radio synchrotron background (RSB) is currently unknown. Its understanding might have profound implications in fundamental physics or might reveal a new class of radio emitters. In this work, we consider the scenario in which the RSB is due to extragalactic radio sources and measure the angular cross-correlation of LOFAR images of the diffuse radio sky with matter tracers at different redshifts, provided by galaxy catalogs and CMB lensing. We compare these measured cross-correlations to those expected for models of RSB sources. We find that low-redshift populations of discrete sources are excluded by the data, while higher redshift explanations are compatible with available observations. We also conclude that at least 20% of the RSB surface brightness level must originate from populations tracing the large-scale distribution of matter in the universe, indicating that at least this fraction of the RSB is of extragalactic origin. Future measurements of the correlation between the RSB and tracers of high-redshift sources will be crucial to constraining the source population of the RSB.

High energy gamma-ray sources in the VVV survey - II. The AGN counterparts

ABSTRACT We identified Active Galactic Nuclei (AGN) candidates as counterparts to unidentified gamma-ray sources (UGS) from the Fermi-LAT Fourth Source Catalogue at lower Galactic latitudes. Our methodology is based on the use of near- and mid-infrared photometric data from the VISTA Variables in the Vía Láctea (VVV) and Wide-field Infrared Survey Explorer (WISE) surveys. The AGN candidates associated with the UGS occupy very different regions from the stars and extragalactic sources in the colour space defined by the VVV and WISE infrared colours. We found 27 near-infrared AGN candidates possibly associated with 14 Fermi-LAT sources using the VVV survey. We also found 2 blazar candidates in the regions of 2 Fermi-LAT sources using WISE data. There is no match between VVV and WISE candidates. We have also examined the Ks light curves of the VVV candidates and applied the fractional variability amplitude (σrms) and the slope of variation in the Ks passband to characterise the near-infrared variability. This analysis shows that more than 85 per cent of the candidates have slopes in the Ks passband >10−4 mag/day and present σrms values consistent with a moderate variability. This is in good agreement with typical results seen from type-1 AGN. The combination of YJHKs colours and Ks variability criteria was useful for AGN selection, including its use in identifying counterparts to Fermi γ-ray sources.

TeV pion bumps in the gamma-ray spectra of flaring blazars

Very high-energy (VHE, $E>100$ GeV) observations of the blazar Mrk 501 with MAGIC in 2014 provided evidence for an unusual narrow spectral feature at about 3 TeV during an extreme X-ray flaring activity. The one-zone synchrotron-self Compton scenario, widely used in blazar broadband spectral modeling, fails to explain the narrow TeV component. Motivated by this rare observation, we propose an alternative model for the production of narrow features in the VHE spectra of flaring blazars. These spectral features may result from the decay of neutral pions ($ bumps) that are in turn produced via interactions of protons (of tens of TeV energy) with energetic photons, whose density increases during hard X-ray flares. We explored the conditions needed for the emergence of narrow $ bumps in VHE blazar spectra during X-ray flares reaching synchrotron energies $ 100$ keV using time-dependent radiative transfer calculations. We focused on high-synchrotron peaked (HSP) blazars, which comprise the majority of VHE-detected extragalactic sources. We find that synchrotron-dominated flares with peak energies $ keV can be ideal periods for the search of $ bumps in the VHE spectra of HSP blazars. The flaring region is optically thin to photopion production, its energy content is dominated by the relativistic proton population, and the inferred jet power is highly super-Eddington. Application of the model to the spectral energy distribution of Mrk 501 on MJD 56857.98 shows that the VHE spectrum of the flare is described well by the sum of a synchrotron self-Compton (SSC) component and a distinct $ bump centered at 3 TeV. Spectral fitting of simulated SSC+$ spectra for the Cherenkov Telescope Array (CTA) show that a $ bump could be detected at a 5sigma significance level with a 30-min exposure. A harder VHE gamma -ray spectrum than the usual SSC prediction or, more occasionally, a distinct narrow bump at VHE energies during hard X-ray flares, can be suggestive of a relativistic hadronic component in blazar jets that otherwise would remain hidden. The production of narrow features or spectral hardenings due to $ decay in the VHE spectra of blazars is testable with the advent of CTA.

Magnetic moments of astrophysical neutrinos

We study the impact of neutrino magnetic moments on astrophysical neutrinos, in particular supernova neutrinos and ultra-high energy neutrinos from extragalactic sources. We show that magnetic moment-induced conversion of Dirac neutrinos from left-handed states into unobservable right-handed singlet states can substantially change the flux and flavour composition of these neutrinos at Earth. Notably, neutrinos from a supernova's neutronisation burst, whose flux can be predicted with 𝒪(10%) accuracy, offer a discovery reach to neutrino magnetic moments  ∼ few × 10-13 μB , up to one order of magnitude below current limits. For high-energy neutrinos from distant sources, for which no robust flux prediction exists, we show how the flavour composition at Earth can be used as a handle to establish the presence of non-negligible magnetic moments, potentially down to few× 10-17 μB if the measurement can be performed on neutrinos from a single source. In both cases, the sensitivity strongly depends on the galactic (intergalactic) magnetic field profiles along the line of sight. Therefore, while a discovery is possible down to very small values of the magnetic moment, the absence of a discovery does not imply an equally strong limit. We also comment on the dependence of our results on the right-handed neutrino mass, paying special attention to the transition from coherent deflection by a classical magnetic field to incoherent scattering on individual scattering targets. Finally, we show that a measurement of Standard Model Dirac neutrino magnetic moments, of order 10-19 μB , could be possible under rather optimistic, but not completely outrageous, assumptions using flavour ratios of high-energy astrophysical neutrinos.

3C 273 host galaxy with Hubble Space Telescope coronagraphy

The close-in regions of bright quasars’ host galaxies have been difficult to image due to the overwhelming light coming from quasars. With coronagraphic observations in visible light using the Space Telescope Imaging Spectrograph (STIS) on the Hubble Space Telescope, we removed 3C 273 quasar light using color-matching reference stars. The observations revealed the host galaxy from 60″ to 0.​​″2 with nearly full angular coverage. Isophote modeling has revealed a new core jet, a core blob, and multiple smaller-scale blobs within 2.​​″5. The blobs could potentially be satellite galaxies or infalling materials towards the central quasar. Using archival STIS data, we constrained the apparent motion of its large scale jets over a 22 yr timeline. By resolving the 3C 273 host galaxy with STIS, our study validates the use of coronagraphs on extragalactic sources for obtaining new insights into the central (at ∼kpc scales) regions of quasar hosts.

Optical spectroscopy of blazars for the Cherenkov Telescope Array – III

Context. Blazars, which include BL Lacs and flat-spectrum radio quasars, represent the brightest persistent extragalactic sources in the high-energy (HE; 10 MeV–100 GeV) and very-high-energy (VHE; E > 100 GeV) γ-ray sky. Due to their almost featureless optical/UV spectra, it is challenging to measure the redshifts of BL Lacs. As a result, about 50% of γ-ray BL Lacs lack a firm measurement of this property, which is fundamental for population studies, indirect estimates of the extragalactic background light, and fundamental physics probes (e.g., searches for Lorentz-invariance violation or axion-like particles). Aims. This paper is the third in a series of papers aimed at determining the redshift of a sample of blazars selected as prime targets for future observations with the next generation, ground-based VHE γ-ray astronomy observatory, Cherenkov Telescope Array Observatory (CTAO). The accurate determination of the redshift of these objects is an important aid in source selection and planning of future CTAO observations. Methods. Promising targets were selected following a sample selection obtained with Monte Carlo simulations of CTAO observations. The selected targets were expected to be detectable with CTAO in observations of 30 h or less. We performed deep spectroscopic observations of 41 of these blazars using the Keck II, Lick, SALT, GTC, and ESO/VLT telescopes. We carefully searched for spectral lines in the spectra and whenever features of the host galaxy were detected, we attempted to model the properties of the host galaxy. The magnitudes of the targets at the time of the observations were also compared to their long-term light curves. Results. Spectra from 24 objects display spectral features or a high signal-to-noise ratio (S/N). From these, 12 spectroscopic redshifts were determined, ranging from 0.2223 to 0.7018. Furthermore, 1 tentative redshift (0.6622) and 2 redshift lower limits at z > 0.6185 and z > 0.6347 were obtained. The other 9 BL Lacs showed featureless spectra, despite the high S/N (≥100) observations. Our comparisons with long-term optical light curves tentatively suggest that redshift measurements are more straightforward during an optical low state of the active galactic nucleus. Overall, we have determined 37 redshifts and 6 spectroscopic lower limits as part of our programme thus far.

Powerful Radio Sources in the Southern Sky. III. First Results of the Optical Spectroscopic Campaign

We recently built the G4Jy-3CRE catalog of extragalactic radio sources. This catalog lists 264 powerful radio sources selected with similar criteria to those of the revised Third Cambridge Catalog, but visible from the Southern Hemisphere. A literature search revealed that 119 sources in the G4Jy-3CRE catalog (i.e., 45%) lack a firm spectroscopic redshift measurement. Here, we present a campaign aimed at acquiring optical spectra of G4Jy-3CRE sources and measuring their redshifts. We used single-slit observations obtained with the Víctor Blanco Telescope, the New Technology Telescope, the Southern Astrophysical Research Telescope, and the 2.1 m telescope of the Observatorio Astronómico Nacional at San Pedro Mártir, Mexico. In addition, we analyzed Very Large Telescope/MUSE archival observations. From these observations, we report the spectra and redshifts of 93 sources, 42 of which are the first optical spectra and redshift determinations for the respective sources. With our new data, approximately 71% of the sources in the G4Jy-3CRE catalog now have firm spectroscopic redshift measurements. This data set will be the basis of our future analysis of the optical properties of the G4Jy-3CRE catalog.

Reduction of ultra-high energy neutrino flux for spacelike neutrinos from extragalactic sources

Assuming that neutrinos are spacelike (tachyonic) fermions, we calculate width for the kinematically allowed, lepton number conserving, three-body decay να→νανβν¯β in the Standard Model. Decays of tachyonic neutrinos over cosmological distances can lead to a reduction of the neutrino flux in the high-energy end of the spectrum. We estimate upper limits on the spacelike neutrino mass based on the PeV-energy cosmological neutrino events observed in the IceCube experiment. These limits are close to those deduced from the measurements of mν2 in the tritium-decay experiment KATRIN.

VLBI position variability of AGNs is inversely correlated with their photometric variability

The stability of the International Celestial Reference Frame (ICRF), realized through geodetic very long baseline interferometry (VLBI) positions of thousands of extragalactic objects, is dependent on the individual positional stability of these objects. It has been recently shown that the prevalence of offsets between the VLBI positions of ICRF objects and their Gaia optical positions, which limit the optical-radio reference frame tie, is inversely correlated with optical photometric variability, suggesting that photometrically variable objects may be more positionally stable. In this work, we determine the relationship between VLBI position stability of ICRF objects and optical-radio position offsets as well as optical photometric variability. We created multi-epoch geodetic VLBI solutions for a sample of 520 ICRF sources that have sufficient data to determine the variability in their VLBI positions over time. We compared this position variability with the fractional photometric variability provided by the Gaia extragalactic source catalog, the Gaia-ICRF optical-radio position offsets, the uncertainty-normalized position offsets, and optical $ BP-RP $ color as well as with possible confounders such as optical magnitude, VLBI/Gaia position error, and redshift. We determined the relationship between VLBI position stability and gamma -ray detection by the Fermi Large Area Telescope (LAT), and we determined how the VLBI position and optical flux variabilities correlate with the spectral classification of our sample, considering flat spectrum radio quasars (FSRQs), quasi-stellar objects, BL Lacs, Seyfert, and gigahertz-peaked spectrum radio sources or compact-steep-spectrum radio sources. We found that VLBI astrometric variability is (i) negatively correlated with optical flux variability, (ii) positively correlated with optical-radio offsets, (iii) negatively correlated with optical color index BP-RP, and (iv) negatively correlated with gamma -ray detection. We also found that the most positionally stable sources are among the FSRQ and BL Lac classes. In other words, redder, photometrically variable sources have the most stable VLBI positions, the smallest optical-radio position offsets, and the highest rate of gamma -ray detection, and these sources tend to be spectrally classified as blazars. Our results are consistent with the most positionally stable sources being blazars, a class of object in which the jet is oriented close to the line of sight and where relativistic beaming increases photometric variability and minimizes the projected offset between the optical and radio positions. Our study should therefore orient future geodetic VLBI observing programs preferentially toward sources with high photometric variability because these sources are predicted to have better VLBI position stabilities and smaller optical-radio position offsets, improving the stability of the celestial reference frame axes.

Influence of sources with a spectral peak in the detection of Cosmic Dawn and Epoch of Reionization

ABSTRACT Foreground removal is one of the biggest challenges in the detection of the Cosmic Dawn (CD) and Epoch of Reionization (EoR). Various foreground subtraction techniques have been developed based on the spectral smoothness of foregrounds. However, the sources with a spectral peak (SP) at megahertz may break down the spectral smoothness at low frequencies (<1000 MHz). In this paper, we cross-match the GaLactic and Extragalactic All-sky Murchison Widefield Array (GLEAM) extragalactic source catalogue with three other radio source catalogues, covering the frequency range 72 MHz–1.4 GHz, to search for sources with spectral turnover. 4423 sources from the GLEAM catalogue are identified as SP sources, representing ≈3.2 per cent of the GLEAM radio source population. We utilize the properties of SP source candidates obtained from real observations to establish simulations and test the impact of SP sources on the extraction of CD/EoR signals. We statistically compare the differences introduced by SP sources in the residuals after removing the foregrounds with three methods, which are polynomial fitting, principal component analysis, and fast independent component analysis. Our results indicate that the presence of SP sources in the foregrounds has a negligible influence on extracting the CD/EoR signal. After foreground subtraction, the contribution from SP sources to the total power in the two-dimensional power spectrum within the EoR window is approximately 3–4 orders of magnitude lower than the CD/EoR signal.

Searching for dark matter annihilation with IceCube and P-ONE

We present a new search for weakly interacting massive particles utilizing ten years of public IceCube data, setting more stringent bounds than previous IceCube analysis on massive dark matter to neutrino annihilation. We also predict the future potential of the new neutrino observatory, P-ONE, showing that it may even exceed the sensitivities of Fermi-LAT gamma-ray searches by about 1–2 orders of magnitude in 1–10 TeV regions. This analysis considers the diffuse dark matter self-annihilation to neutrinos via direct and indirect channels, from the galactic dark matter halo and extra-galactic sources. We also predict that P-ONE will be capable of pushing these bounds further than IceCube, even reaching the thermal relic abundance utilizing a galactic center search for extended run-time.

Constraints on the proton fraction of cosmic rays at the highest energies and the consequences for cosmogenic neutrinos and photons

Over the last decade, observations have shown that the mean mass of ultra-high-energy cosmic rays (UHECRs) increases progressively toward the highest energies. However, the precise composition is still unknown and several theoretical studies hint at the existence of a subdominant proton component up to the highest energies. Motivated by the exciting prospect of performing charged-particle astronomy with ultra-high-energy (UHE) protons we quantify the level of UHE-proton flux that is compatible with present multimessenger observations and the associated fluxes of neutral messengers produced in the interactions of the protons. We study this scenario with numerical simulations of two independent populations of extragalactic sources and perform a fit to the combined UHECR energy spectrum and composition observables, constrained by diffuse gamma-ray and neutrino observations. We find that up to of order 10% of the cosmic rays at the highest energies can be UHE protons, although the result depends critically on the selected hadronic interaction model for the air showers. Depending on the maximum proton energy (E max p) and the redshift evolution of sources, the associated flux of cosmogenic neutrinos and UHE gamma rays can significantly exceed the multimessenger signal of the mixed-mass cosmic rays. Moreover, if E max p is above the GZK limit, we predict a large flux of UHE neutrinos above EeV energies that is absent in alternate scenarios for the origin of UHECRs. We present the implications and opportunities afforded by these UHE proton, neutrino and photon fluxes for future multimessenger observations.

Lost in the curve: Investigating the disappearing knots in blazar 3C 454.3

One of the most well-known extragalactic sources in the sky, quasar 3C 454.3, shows a curved parsec-scale jet that has been exhaustively monitored with very-long-baseline interferometry (VLBI) over the recent years. In this work, we present a comprehensive analysis of four years of high-frequency VLBI observations at 43 GHz and 86 GHz, between 2013–2017, in total intensity and linear polarization. The images obtained from these observations enabled us to study the jet structure and the magnetic field topology of the source on spatial scales down to 4.6 parsec in projected distance. The kinematic analysis reveals the abrupt vanishing of at least four new superluminal jet features in a characteristic jet region (i.e., region C), which is located at an approximate distance of 0.6 milliarcsec from the VLBI core. Our results support a model in which the jet bends, directing the relativistic plasma flow almost perfectly toward our line of sight, co-spatially with the region where components appear to stop.

Enlightening the Universe behind the Milky Way bulge

Context. The location of the Solar System complicates the detection extragalactic sources beyond the Milky Way plane. The optical observations are hampered in the so-called zone of avoidance (ZOA), where stellar crowding and Galactic absorption are severe. Observations at longer wavelengths are needed in order to discover new background galaxies and complete the census in the ZOA. Aims. The goal of this work is to identify galaxies behind the Milky Way bulge using near-infrared (NIR) data from the VISTA Variables in Vía Láctea (VVV) survey. Methods. To this end, we made use of different VISTA Science Archive (VSA) tools in order to extract relevant information from more than 32 billion catalogued sources in the VVV bulge region. We find that initial photometric restriction on sources from the VSA vvvSource table combined with restrictions on star–galaxy separation parameters obtained from Source Extractor is a successful strategy for achieving acceptable levels of contamination (60%) and high completeness (75%) in the construction of a galaxy target sample. To remove contaminating Galactic sources from the initial target sample, our methodology also incorporates a visual inspection of false-colour RGB images, a crucial quality control carried out following a specifically defined procedure. Results. Under this methodology, we find 14 480 galaxy candidates in the VVV bulge region, making this sample the largest catalogue to date in the ZOA. Moreover, these new sources provide a fresh picture of the Universe hidden behind the curtain of stars, dust, and gas in the unexplored Milky Way bulge region. Conclusions. The results from this work further demonstrate the potential of the VVV/VVVX survey to find and study a large number of galaxies and extragalactic structures obscured by the Milky Way, expanding our knowledge of the Universe in this challenging and impressive region of the sky.